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Abstract 

We announce a detailed numerical investigation for some class of difference schemes, which arises 
from Euler implicit scheme. Such schemes demonstrate unusual behavior and leads to origin of 
solitons. Applications to some nonlinear problems are discussed. 
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For the numerical calculations in mathematical physics, one usually applies standard 
procedure in which initial differential equation is replaced by system of finite-difference 
equations in accordance with some difference scheme. However, contrary situations when 
the problem is formulates initially in the terms of some difference scheme also are known. 
May be the best example is the numerical FPU experiment. In 1955, Fermi, Pasta, and 
Ulam took numerical scheme for the chain of harmonic oscillators coupled with a quadratic 
nonlinearity and investigated how the energy in one mode would spread to the rest 
They found the system cycled periodically, implying it was much more integrable than 
they had thought. Only in 1963, Martin Kruskal and Norman Zabusky analytically study 
the Fermi-Pasta-Ulam heat conduction problem in the continuum limit and find that the 
Korteweg - de Vries equation governs this system Q]. 

In our note, we consider similar situation when the problem with very intriguing solu- 
tions, which was formulates for the some different scheme, has not continual analogue yet. 
Discussed problem comes from the rounding error analysis of the Euler implicit scheme for 
the linear diffusion equation 

Ui,j+1 (1 + 2 ■ c) + C • Mi-lj+l + C ■ Ui+ij+l = u id . 

Let's replace in this scheme 2 with some arbitrary coefficient A 

Uij+i (1 + A • c) + c ■ Uj-ij+i + c • = Uij (1) 

and examine numerically how solutions will changes with changing of A. Let's take initial 
conditions as ua t i = 1 and boundary conditions as u\j = and u n j = 0, where % G [l,n] 
and j G [1, m]. 

Obtained results are as follows. For A = 2 and for real c one obtain standard solution for 
the diffusion equation (Fig. [IJ, for the complex c - usual solution of the linear Schrodinger 
equation (Fig. EJ). 

For the A > 2 solution stay stable for both real and imaginary values of c (Fig. |3] and 
Fig. IU). 

However, results for A < 2 seems to be very different. For all real c the family of schemes 
(1) gives divergent solutions (jHJ). 

More intriguing solutions arises in the case of imaginary c. Namely, in this situation for 
all real values from interval A G (0,2) soliton-like formations appear. Some results for this 
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FIG. 1: Solution for the linear diffusion equation, c = —1, i G [1,500], j € [1,30] 



case we present in Fig. |B] Fig. |H| In Fig. H where A = 1.9 for the closest "time layers" j 
solitons has bad resolution. 

At the same time solitons velocities and its sharpness grows up with decreasing of A. In 
Fig. 0we present results of calculations for A = 1 and the same value of c. 

As a whole soliton's behavior is similar to those of solitons of the Korteweg - de Vries 
equation. Anyway, its form, character of reflection from the limit points of the integration 
interval, its scattering all those factors are very and very similar to KdV solitons (Fig. |H] 
and FigEI). With further decreasing of A soliton's velocity increased and one can observe its 
evolution for the deep "time" layers j. 

In the discussed family of difference schemes solitons occurs only in the case when the 
initial conditions contains the discontinuity of u. At that, amplitude is proportional to 
gradient of u but velocity of solitons stays the same. One can see those features in Fig. [TU] 
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and u — 1 for i £ 



+ l,n 



and FiglTTl where initial conditions is u = 0.9 for i 6 

Presented results of calculations for A < 2 seems to be amusingly enough but it is a 
question how those schemes and anomalies correlates with known mathematical (and may 
be physical) objects. Probably, it is the task of further investigations but some speculations 
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FIG. 2: Solution for the linear Schrodinger equation, c = —i 
are well-timed now. 

Firstly, one can note that soliton's origin realizes in anyone difference scheme where A < 2. 
If A weakly differ from 2 effects occurs in very deep "time" layers (for the big j-s) but we 
believe that effect actual exists. Therefore, such a circumstance one may took in account 
even for the standard scheme with A = 2. 

Next, it is incomprehensible now what kind of problems can correspond to schemes with 
A < 2 and only rough assumptions are permissible. Though founded solutions similar in 
appearance with KdV solitons they appear only for imaginary values of c. Therefore, one 
can assume that such solutions relates to Schrodinger equation with one or another kind of 
nonlinearity. It seems to be very stirring. 

Indeed, let's consider the following problem. Let A = 1 and initial conditions are given 
in the form of some harmonic function as it shown in FigEl • Let narrow barrier divides 
left and right parts of the integration interval. Let's put initial function inside the barrier 
and on the right of barrier identical with zero. 

Now, let's observe evolution of the initial function with the "time". In the left border 
of interval and in the left boundary of barrier solitons appear. Soliton 1 became its motion 
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FIG. 3: Solution for the difference scheme with parameters c = — 1, A = 2.1 

to the right, soliton 2 to the left, and soliton 3 to the right. Solitons 1 and 2 collides 
with each other and soliton 3 continues its motion. However, it is notable that ratio of 
integrals of squares of absolute values of considered function stays exactly invariable for the 
long "time". In the case of A = 2 such a behavior should corresponds to the underbarrier 
tunneling and above ratio one may consider as tunneling rate (see Fig|13|). However, in our 
case it is nothing more than interesting peculiarity. 

We hope that further investigations in this direction may be perspective enough. 
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FIG. 7: Solitons for c = — j, A = 1 



FIG. 8: Scattering of solitons, c = —i, A = 1, m = 350 
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FIG. 11: Velocities of solitons, which occurs in the limits of integration range and in its center, are 
equal 
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FIG. 12: "Soliton-type" tunneling 
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FIG. 13: "Tunneling" rate for A = 1. Origin of the plateau corresponds with escape of soliton 
3 from under the barrier; ending of the plateau corresponds with entering of soliton 3 under the 
barrier 
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